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Using molecular dynamics simulations, we show that, when subject to a periodic external elec- 
tric field, a nanopore in ionic solution acts as a capacitor with memory (memcapacitor) at various 
frequencies and strengths of the electric field. Most importantly, the hysteresis loop of this mem- 
capacitor shows both negative and diverging capacitance as a function of the voltage. The origin 
of this effect stems from the slow polarizability of the ionic solution due to the finite mobility of 
ions in water. We develop a microscopic quantitative model which captures the main features we 
observe in the simulations and suggest experimental tests of our predictions. We also suggest a 
possible memory mechanism due to the transport of ions through the nanopore itself, which may be 
observed at small frequencies. These effects may be important in both DNA sequencing proposals 
using nanopores and possibly in the dynamics of action potentials in neurons. 



A nanopore is an aperture of nanoscale dimensions 
across an insulating membrane. This way, if immersed in 
an ionic solution, ions can enter the nanopore, and if sub- 
ject to an electric field of given strength, they may enter 
it from one side of the opening and emerge on the other 
side. The membrane can be made of either biological or 
solid-state materials. Nanopores are ubiquitous in biolog- 
ical systems as they regulate the flow of ions across cell 
membranes (of, e.g., neurons). At present, they are also 
actively investigated for potential applications in DNA 
sequencing fT, ^2] . Despite its apparent simplicity, ion dy- 
namics in nanopores is far from trivial. An example of 
this is provided by the recently predicted phenomenon 
of "quantized ionic conductance", namely the presence of 
current steps occurring at effective pore radii that match 
the radii of the water hydration layers that form around 
the ions [3j. 

Here, we predict another intriguing property of ion dy- 
namics in nanopores, when the latter are subject to a 
periodic external electric field. In particular, we show 
that the electric field forces ions to accumulate at the 
two surfaces of the nanopore membrane, creating an ef- 
fective capacitor. However, this capacitor shows inter- 
esting features as a function of the frequency of the field, 
namely a hysteresis loop of the capacitance (memory- 
capacitance or memcapacitance for short [4J) as a func- 
tion of voltage which diverges at zero voltage and dis- 
plays negative-capacitance properties. Memcapacitors 
have been theoretically formulated in Ref. [4j and be- 
long to the wider class of memory-circuit elements, which 
includes also memristors and meminductors, namely re- 
sistors and inductors with memory. Following Ref. [4J we 
define a voltage-controlled memcapacitive system by the 
set of equations 

Q{t) = C{x,V,t)V{t) (1) 
x = f{x,V,t) (2) 



where Q{t) is the charge on the capacitor at time t, V{t) is 
the voltage across it, and C is the time-dependent mem- 
capacitance which depends on the internal state of the 
system described by a set of n state variables x, with / 
a continuous n-dimensional vector function. Using all- 
atom molecular dynamics (MD) simulations, we indeed 
demonstrate that memory characteristics of nanopores 
in solution originate from the finite mobility of ions in 
water with consequent slow polarizability of the ionic so- 
lution. We also develop a simple microscopic model that 
captures the main properties observed in the simulations 
and allows us to extend our results to regimes beyond 
the reach of MD simulations. Additionally, we use an 
equivalent circuit formulation to discuss how ionic trans- 
port through the nanopore itself may lead to an addi- 
tional memory mechanism. We also propose ways our 
predictions could be tested with available experimental 
capabilities. 

Ionic memcapacitors — In order to show that a 
nanopore in solution indeed acts as a memcapacitive sys- 
tem (Eqs. ([T]) and Q), let us be more specific and con- 
sider the experimental situation as depicted in Fig. [l] 
which we suggest as a way to observe the effects pre- 
dicted here. An external (time-dependent) voltage V{t) 
is applied to electrodes which form an effective capaci- 
tor with the ionic solution and the nanopore at its inte- 
rior. In this paper we consider a typical setup for, e.g., 
DNA sequencing devices [SHZl with KCl ionic solution 
and a Si3N4 nanopore 25 A thick and 88 A wide with a 7 
A radius cylindrical hole drilled through the center of the 
membrane. However, the specific dimensions of the pore 
are irrelevant for the overall conclusions of this paper 
and similar considerations would apply also to biologi- 
cal pores. When a global time-dependent bias is applied 
to the system, the ionic solution is polarized and ions 
are forced to accumulate at the surfaces of the pore (see 
Fig. fl]). First, we consider relatively short time scales 
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when ion transport across the nanopore is not signifi- 
cant. (We will discuss towards the end of the paper the 
effect of ion transport across the pore at longer times.) 
This implies that the same (and opposite in sign) charge, 
Q{t), that accumulates on the surfaces of the pore also 
accumulates on the plates of the capacitor. 

In the simplest approximation, the electrical proper- 
ties of the total nanopore system can be approximated 
by an equivalent circuit model as shown in Fig. [l] where 
capacitors and resistors (generally, non-linear) represent 
different parts of the system. Here, Ce-s denotes the ex- 
ternal electrode-solution capacitance, Rs is the resistance 
of the solution, Rgg (^ Rs) is the resistance of the ion 
current through the pore, and Css is the capacitance of 
the membrane. The total voltage drop is then given by 



2R, 



dQ Q_ 

' dt ^ 



V{t). 



(3) 



Since we are interested in the properties of the pore and 
solution only, we envision the experimental set-up so that 
Ce-s ^ Css' The equation for the charge then takes the 
form 



dQ ^ V(t)_ 
dt ~ 2Rs 



2RsCs 



(4) 



Eq. Q describes a relaxation of Q towards V{t)/{2Rs) 
with a relaxation time 2RsCss- Generally, due to the 
finite mobility of ions in water, the charge on the ex- 
ternal plates - and hence on the pore - is slow to re- 
spond to changes of the bias upon varying V{t). This 
means that V{t) may be zero when Q{t) is finite and vice 
versa. Therefore, the capacitance of this whole system, 
C = Q{t)/V{t)^ depends on the applied voltage history. 
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Figure 1. (Color online) Left: a snapshot of the molecular dy- 
namics geometry at a time when a buildup of charges of the 
opposite sign on each side of the nanopore is observed due to 
a finite electric field E. The nanopore membrane is located at 
the center (represented by the red lines) and water is not visi- 
ble. Top and bottom horizontal red lines represent electrodes 
(holding plate charges di(5(t)) of the suggested experimental 
set-up needed to observe the predicted memcapacitive effects. 
Right: simplified equivalent circuit model. 
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Figure 2. (Color online) Net charge on the surface of the 
nanopore (equal to system's plate charge) and voltage across 
the nanopore (which is very close to the external voltage, 
solid line) plotted as a function of time. The solid line for 
the charge corresponds to the model we discuss in the text. 
One can clearly see that the net charge lags behind in time 
with respect to both the external voltage and the voltage at 
the pore. 



shows divergences and acquires negative values at spe- 
cific times, similar to what has been recently predicted 
in some solid-state memcapacitors f8]. Since these mem- 
ory effects can be viewed as originating from the history- 
dependent permittivity of the ionic solution, following 
the definition of Ref. [4J (Eqs. ([T]) and ([2])), a nanopore 
in ionic solution is indeed a memcapacitive system. 

Results and discussion — Let us now demonstrate 
these memcapacitive effects from a microscopic point of 
view. For this we employ all-atom MD simulations using 
NAMD2 [9] and investigate the response of the system 
to external ac- and dc-electric fields. According to our 
previous discussion, in order to calculate the total capac- 
itance of the system we only need the net charge that 
accumulates at the nanopore surface. We then employ 
periodic boundary conditions in both the direction per- 
pendicular and parallel to the external field. On each 
side of the pore, we place a 50 A long section of water 
with a 1 M solution of a homogenous, random distribu- 
tion of potassium and chlorine ions. The CHARMM27 
force field p^Of [TT] is used for the interaction of water and 
ions while quantum mechanical parameters [12] are used 
for the Si3N4 pore [13j. 

We have performed extensive MD simulations of the 
process of polarization of the ionic solution in proximity 
to the nanopore membrane. As a general observation, 
within the parameters used, a delay between the forma- 
tion of ionic polarization and both the external ac-electric 
field and the voltage across the pore has been observed. 
This is evident in Fig. [2] where the net charge of ions on 
one side of the nanopore is plotted as function of time to- 
gether with the voltage drop across the membrane. The 
latter is obtained by doing a full electrostatics calculation 
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from MD simulations [T^. At this point we remark that 
in ah our simulations we have found that the net charge 
induced by the external field is always located within 10 
A from the membrane surface and that the voltage drop 
across the membrane is very close in value to the external 
applied voltage, whatever the strength and frequency of 
the latter. This last observation is consistent with previ- 
ous studies [T]. 

Next, we fit these MD simulation results employing 
the equivalent circuit model given by Eq. Q. Using 
microscopic parameters, Eq. ^ can be rewritten as 



dQ{t) 
dt 



2eAiinE(t) 



Q{t) 



(5) 



where A is the area of the membrane surface, /i is the ion 
mobility which is similar for both types of ions and we 
take to be 7.12 x 10"^ m^l/"^s"^ [15j, e the ion charge, 
n the density of ions in the bulk, and r = 2RsCss is the 
relaxation time. The factor of 2 takes into account the 
conductivity of both channels (K+ and Cl~). Using dc- 
field simulations, we were able to extract r. We do this 
by first applying a constant electric field (with the same 
magnitude as the amplitude of our ac-field simulations) 
which forces ions to accumulate in the vicinity of the 
membrane as it is shown in Fig. [T] The total net charge 
within 10 A from the nanopore surface is integrated to 
obtain the capacitor charge [16j. Then, we turn off the 
field and monitor the decay of this charge back into the 
bulk (see Fig. [3|. From this we obtain the relaxation 
time r. We find that this time is slightly dependent on 
the amplitude of the initially applied electric field, be- 
coming slightly longer as the applied field decreases. For 
voltages as low as 5 V across the membrane we get from 
our simulations a relaxation time r 75ps [17J. 

Solving Eq. ([5| for Q{t) with an ac electric field E{t) = 
Eq sin {uot) of amplitude Eq and angular frequency uj = 
27r/, we obtain the long-time limit solution: 



Q{t) *" 



2AfienEQ sin (ut) 2AfienEQ cos {ut) 



(1 



(1 



(6) 



which amounts to a sine function with a phase shift, 
namely Q{t) = Qo sin {ut — (5), where 
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2AiienEQ 



(1 



1/2 



(7) 

and 5 = tan ^ {uot). Using the appropriate values of 
A, n, and r as obtained from the numerical simulations, 
this model agrees very well with the results of our MD 
simulations as it is evident from Fig. [2] In particular, it 
is clear that the net charge and voltage on the capacitor 
are phase shifted by the amount 5. This means that 
when the voltage is zero the charge on the capacitor is 
not necessarily zero, and vice versa. This is shown in 
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Figure 3. (Color online) Net charge of the positively charged 
side of the the capacitor vs. time when a constant 20 V ap- 
plied voltage responsible for the accumulation of the charges 
on the surface of the pore is turned off. A decay timescale 
of 75 ps is obtained by fitting the simulation results with 
an exponential decay curve, i.e., Q(t) — Qmax^~^^^ - Due 
to the nearly equal mobilities of and Cl~ [15 , we find 
similar results for the negatively charged side and different 
field strengths as well. This is shown in the inset with initial 
voltage of 5 V. 



Fig. [IJa) where the net charge is plotted as a function 
of the voltage across the capacitor. This gives rise to 
diverging and negative values of capacitance. 

To see this, we write the external voltage across the 
whole system as V(i) = Fosin(ct;t), with = —E{^d 
where d is the size of the simulation cell along the direc- 
tion of the electric field. From C = Q{t)/V(t) and ^ 
we then obtain the long-time limit of the capacitance: 



C(t) 



2Ajiien 



2Ajj.en 



(1 



j^Td (1 



jd 



cot {cot), 



(8) 

This represents the main result of our paper. It shows 
that the capacitance of the whole nanopore system is 
history-dependent, diverges when the voltage crosses its 
zero value, and acquires negative values within a certain 
voltage range. This is plotted in Fig. [4] where both the 
charge on the capacitor and the capacitance are plot- 
ted as a function of voltage across the capacitor. From 
Fig. [4] it is evident that at low frequencies, this ionic 
memcapacitor behaves almost like a regular capacitor. 
At higher frequencies, memory effects in the capacitance 
manifest in a hysteresis loop characteristic of memcapac- 
itors. However, unlike typical memcapacitors that show 
a pinched hysteretic loop [4J, ionic memcapacitors have 
a non- vanishing (diverging) zero-bias capacitance. 

Memory at small frequencies — At small enough fre- 
quencies, there will be enough ion transport through the 
pore to play a role in the dynamics. In fact, it is no 
longer correct to assume that the charge associated with 
Ce-s is the same as the charge associated with Css- Us- 
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Figure 5. (Color Online) Memcapacitive effects due to ionic 
transport across the pore. Main panel shows the charge on 
the external capacitor plates plotted versus the voltage across 
it for various frequencies of the electric field with Vb = 1-0 V, 
A = 10 /xm^, Lp = 25 A, rp = 7 A, e = 7.5eo, Css eA/Lp, 
and Ce-s = lOCss for a 1 M concentration of ions. The inset 
shows the capacitance plotted versus the voltage for the same 
frequencies. Memory effects due to ionic transport across the 
pore occur at much lower frequencies than those due to the 
polarization of the ionic solution. 




Figure 4. (Color online) Upper panel: net charge Q versus 
a periodic voltage of amplitude Vb = 1 V and different fre- 
quencies as obtained from Eq. |6] The inset shows the same 
quantity for 2 GHz and Vb = 20 V compared with the charge 
obtained directly from our MD simulations. Lower panel: ca- 
pacitance versus a periodic voltage of frequency / =10 MHz 
and amplitude Vb = 1 V as obtained from Eq. |8] In the in- 
set we show the same quantity for / =2 GHz and Vb = 20 
V compared with the capacitance obtained directly from our 
MD simulations. It is seen that C can be both negative and 
diverges as the voltage approaches zero. The arrows indicate 
the direction the voltage is swept in time and the numbers 
show the order in which the trace is generated. 



ing Kirchoff's loop rule, we get an equation similar to 
Eq. (§: 



2R 



dQe 



dt 



V{t). 



(9) 



except now we differentiate between the two charges, 
Qe-s {on the external plates) and Qss (on the pore). 
Applying Kirchoff's junction rule we then get 



dQe 



dt 



dQs, 
dt 



(10) 



where due to Css and Ess being in parallel, Vss is sim- 
ply the voltage on the Css capacitor, Qss/Css^ and we 
can approximate the conductance across the pore as 
l/Rss = 2Apenfi/Lp where Ap = 7rr^, with Tp the pore 



radius, e the electronic charge, n the ion density, /i the 
ion mobility, and Lp the length of the nanopore. We ap- 
proximate Css as that of a parallel plate capacitor, i.e., 
Css = ^A/Lp with e the permittivity of Si3N4, and A is 
again the area of the membrane surface. We can solve 
Eq. [9] for 



and plug it into Eq. 
following second order equation for Q 



10 to arrive at the 
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dt 



' a,^ (11) 



R 



dt 



For an ac field, this equation can be solved analytically 
but the solution is a bit involved. We then just plot it for 
different frequencies in Fig. [5] As anticipated, we see that 
direct transport across the pore leads to an additional 
memory mechanism at specific frequencies, much smaller 
than the ones due to the polarization of the ionic solution. 

Experimental test — The results presented in this pa- 
per show that there are two very distinct memory regimes 
for nanopores in solution. To observe the first memory 
mechanism due to the polarization of the ionic solution, 
ideally, one would need high frequencies (on the order 
of IGHz). On the other hand, much smaller frequencies 
(on the order of IKHz) are necessary to observe the addi- 
tional memory effect due to transport through the pore. 
This indicates that the second memory mechanism would 
be the easiest one to detect. However, in both cases we 
have shown that the capacitance acquires diverging (and 
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also negative) values. This result appears to hold also 
at frequencies that are relatively smaller than those re- 
quired to observe the first memory regime, namely at 
frequencies that are on the order of tens of MHz or less 
(see Fig. [4|. Therefore, with appropriate control of the 
external circuit, at these frequencies one should at least 
be able to observe non-trivial changes (manifested, e.g., 
in fast jumps between positive and negative values) of 
the capacitance when the bias crosses its zero value. We 
thus expect that by simultaneously measuring the elec- 
tric charge on the external electrodes as shown in Fig. [T] 
(or, equivalently, the current at the electrodes) in the 
presence of an applied time-dependent voltage should al- 
low a direct verification of our predictions for both types 
of memory. 



Conclusions — We have shown, using molecular dy- 
namics simulations, that a nanopore sequencing setup 
acts as a memcapacitor, namely a capacitor with mem- 
ory j4j. The latter is due to two types of effects, and thus 
arises at very different frequencies of the external bias. 
At high frequencies the finite mobility of ions in water 
and hence the slow polarizability of the ionic solution 
give rise to one type of memory. Memcapacitive effects, 
however, may also occur as a result of the ion transport 
through the pore at very low frequencies. These pro- 
cesses occur internally in the system and, from the point 
of view of an external circuit, the whole system behaves 
as an unusual capacitor. These effects may potentially 
play a role in nanopore DNA sequencing proposals, es- 
pecially those based on ac-electric fields [18], as well as 
in other nanopore sensing applications. Moreover, the 
effect of the charge buildup on the nanopore surface may 
influence DNA translocation and its structure in proxim- 
ity to the pore. Finally, due to the ubiquitous nature of 
nanopores in biological processes, these results may be 
relevant to specific ion dynamics when time-dependent 
fields are of importance, such as in the action potential 
formation and propagation during neuronal activity. We 
thus hope this work will motivate studies in this direc- 
tion. 
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